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This  investigation  was  performed  by  Sybucon,  Inc.,  under 
Contract  No.  DAAG29-76-C-0045,  for  the  U.  S.  Army  Research 
Office,  Research  Triangle  Park,  N.  C.  27709.  The  Technical 
Monitor  was  Dr.  Robert  E.  Singleton.  Dr.  John  F.  Nash  and 
Dr.  Roy  M.  Scruggs  served  as  joint  Principal  Investigators. 


ABSTRACT 


Fundamental  investigations  have  been  performed  on  the  behavior  of 
unsteady  compressible  turbulent  and  laminar  boundary- layer,  and 
boundary-layer-like,  flows  experiencing  local  reversal  and  separation. 
An  advanced  numerical  computational  method  was  used,  based  on  an 
implicit-ADI  numerical  integration  technique.  The  computing  capa¬ 
bility  and  economy  of  the  method  permit  a  systematic  study  of  the 
subject  flows,  from  first-order  to  higher-order  representations.  The 
research  extended  earlier  work  performed  by  Sybucon  which  was  also 
supported  by  the  Army  in  collaboration  with  NASA  Ames  Research  Center. 


NOTATION 


h 

P 

t 

u,v,w 

x,y,z 

Aj  through  A^ 

F 

G 

H 

L 

P.Q 

Q 

T 

u.v.w 

p 

x’  z 
e 

V 

Bar  O 


total  heat  flux 

pressure 

time 

fluctuating  velocity  components 
local  orthogonal  coordinates 
numerical  integration  matrices 
vector  of  dependent  variables 
vector  of  values  of  F 
total  enthalpy 

turbulence  dissipation  length 
numerical  integration  matrices 
total  mean  velocity 
temperature 

mean  velocity  components 
density 
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I  INTRODUCTION 


Unsteady  boundary- layer  separation  is  a  limiting  factor  in  the  performance 
of  many  aerodynamic  systems:  helicopter  rotors,  turbine  blades,  the  lifting- 
surfaces  of  high-maneuverability  vehicles,  and  so  forth  [1].  The  successful 
design  of  such  systems  is  far  from  straightforward  because  of  the  complexity 
of  the  associated  fluid  mechanics,  and  because  experimental  data  are  scarce 
and  costly  to  obtain.  In  the  absence  of  comprehensive  measurements,  know¬ 
ledge  of  unsteady  separating  flows  is  being  acquired  largely  through  the 
intelligent  use  of  numerical  experimentation.  Again  because  of  the  lack  of 
comprehensive  measurements,  the  numerical  experiments  are  being  performed 
using  models  whose  validity  can  only  be  inferred  by  extrapolation  of  results 
from  other  flow  regimes:  primarily  steady  flows.  The  need  for  this  extra¬ 
polation  places  increased  emphasis  on  the  need  for  fundamental  correctness 
in  the  models  themselves,  and  a  certain  sophistication  of  the  models  is 
called  for,  subject  to  the  obvious  economic  constraints. 

Extensive  calculations  have  been  carried  out,  over  the  last  several  years,  to 
investigate  the  properties  of  time-dependent  turbulent  boundary  layers  [2,3,4, 
5,6].  It  has  been  shown  that  the  effects  of  time-dependence  are  to  delay  the 
onset  of  flow  reversal  in  the  boundary  layer  [4],  and,  in  many  cases,  to 
delay  the  onset  of  separation  --  in  the  sense  of  detachment  of  the  outer  flow 
from  the  body  surface  --  even  more  [5,6].  Reversal  and  separation  are  dis¬ 
tinct  events,  in  unsteady  flow  [7],  and  between  them  lies  a  flow  regime  in 
which  the  boundary  layer  remains  thin  despite  the  presence  of  reversed  flow 
close  to  the  body  surface.  The  existence  of  this  intermediate  regime  has  an 


important  bearing  on  the  sequence  of  events  which  lead  up  to  dynamic  stall 
[1,8],  and  it  permits  the  use  of  first-order  boundary- layer  theory  to  study 
the  flow  characteristics  therein  [5,6].  Studies  of  this  latter  type  have 
been  made  for  flows  which  approach  separation  monotonically  [4,5],  for  flows 
which  are  perturbed  but  subsequently  allowed  to  relax  towards  steady-state 
conditions  [9],  and  for  oscillatory  flows  in  which  the  separation-provoking 
conditions  are  alleviated  during  part  of  each  cycle  [9].  In  flows  subject 
to  a  progressively  severe  retardation  of  the  outer  inviscid  stream,  a  region 
of  reversal  forms  and  grows  in  extent  with  increasing  time.  As  the  reversal 
develops,  the  displacement  thickness  grows  larger,  but  it  is  not  until  later 
that  a  more  rapid  increase  of  displacement  thickness  accompanies  the  onset 
of  a  singularity  in  the  solution.  Interesting  and  important  data  were 
generated  by  all  of  these  studies,  leading  to  an  increasing  body  of  know¬ 
ledge  about  the  reversal/separation  phenomenon. 

With  the  intention  of  elucidating  the  details  of  the  flow  in  the  vicinity 
of  separation,  some  calculations  were  done  for  unsteady  laminar  flow  using 
the  two-dimensional,  incompressible  Navier-Stokes  equations  [9],  The  Navier- 
Stokes  equations  do  not  become  singular,  when  separations  occurs,  but  do 
portray  the  increasing  displacement  thickness  in  situations  where  a 
singularity  appears  in  the  first-order  boundary- layer  equations.  The 
results  of  the  calculations  using  the  Navier-Stokes  equations  helped  to 
clarify  the  solutions  of  the  turbulent  boundary- layer  equations  obtained 
earlier,  and  also  served  to  relate  the  work  on  turbulent  flows  at  Sybucon 
to  studies  made  elsewhere  for  purely  laminar  flows  [10,  eg.]. 
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The  work  done  by  Sybucon  on  time-dependent  turbulent  boundary  layers,  up  to 
the  time  of  the  present  study,  was  all  done  using  a  calculation  method  based 
on  the  integration  of  the  partial  differential  equations,  together  with  a 
turbulence  model  involving  the  turbulent  kinetic-energy  equation,  by  means 
of  an  explicit  numerical  scheme.  Although  the  method  served  to  provide  a 
wide  range  of  useful  results,  significant  weaknesses  became  apparent.  One 
of  these  was  the  necessity  for  matching  the  numerical  solution  to  the  Law 
of  the  Wall  in  the  region  close  to  the  body  surface.  The  use  of  approximate 
relationships  in  the  wall  region  --  which  is  crucially  important  to  the 
dynamics  of  the  reversed  flow  --  was  of  increasing  concern.  Another  weak¬ 
ness  was  the  use  of  an  explicit  numerical  scheme  which  leads  to  excessively 
long  computer  run-times  for  flows  of  practical  interest.  Furthermore,  the 
internal  determination  of  integration  step  size,  necessitated  by  the 
explicit  formulation,  made  it  more  difficult  to  describe  certain  types  of 
flow  --  notably  oscillatory  flows  --  than  would  be  the  case  if  the  temporal 
node  points  could  be  specified  freely  by  the  user.  Another  limitation  in  the 
existing  method  (not,  this  time,  related  to  the  explicit  numerical  scheme) 
was  its  restriction  to  incompressible  flows.  Aerodynamic  systems  rarely 
operate  in  purely  incompressible  flow,  and  it  was  felt  to  be  of  interest  to 
explore  the  effects  of  compressibility  on  unsteady  reversal  and  separation. 
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II  THE  CALCULATION  METHOD 


The  method  involves  the  integration  of  the  following  governing  equations, 
in  which  x  and  2  are  measured  in  the  developed  surface  of  the  cylinder  (x 
normal  to  the  generators,  and  2  along  them),  and  y  is  measured  normal  to 
the  developed  surface: 
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together  with  the  continuity  equation: 


r—  (pU)  +  (pV)  =  0 
Dx  dy 


(5) 


the  definition  of  the  convective  derivative: 
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and  the  usual  peripheral  relationships: 
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The  turbulent  shear-stress  components  are  assumed  to  be  functions  of  the 
turbulent  kinetic  energy  and  the  mean  rate  of  strain: 


uv 
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where  L  is  the  dissipation  length,  and  the  corresponding  relationship  for 
the  turbulent  heat  flux  is 


V(j) 


3T 

By 


(16) 
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The  dissipation  length  is  assumed  to  be  the  usual  prescribed  function  of 
y/6£;  however,  is  not  now  the  local  boundary- layer  thickness  itself,  but 
is  a  length  related  to  it  via  a  rate  equation  whose  characteristic  time  is 
66/Qg.  The  effect  of  the  introduction  of  this  rate  equation  is  to  place  the 
turbulence  model  on  essentially  the  same  status  as  the  two-equation  models 
used  by  several  other  investigators. 

Experimental  data  for  steady  flow  indicates  that  the  function  f,  in  Equations 
(14,15),  can  be  assumed  to  be  linear.  The  corresponding  form  for  g,  in 
Equation  (16),  can  be  determined  from  f  by  specifying  the  distribution  of 
turbulent  Prandtl  number  across  the  boundary  layer.  In  the  present  cal¬ 
culations  this  Prandtl  number  is  taken  to  be  a  constant,  but  the  method  is 
not  restricted  to  this  form. 

The  above  equations  are  formulated  for  turbulent  flow.  Calculations  for 
laminar  flow  are  performed  simply  by  suppressing  the  turbulent  kinetic 
energy,  i.e.,  by  specifying 

~2 

q  =0, 

whereupon  the  equations  become  identical  to  the  compressible  laminar  boundary- 
layer  equations.  It  will  be  noted  that  the  viscous  terms  are  retained  in  the 
momentum  and  thermal  energy  equations.  The  solution  for  turbulent  flow  is 
continued  through  the  viscous  sublayer  to  the  body  surface;  the  solution  is 
not  matched  to  a  separate  inner-layer  calculation  in  the  wall  region. 

The  four  principal  governing  equations  (2  momentum,  thermal  energy,  and 


turbulent  kinetic  energy)  can  conveniently  be  written  as  a  matrix-vector 
equation: 


A1F  + 


A  + 
2  3t 


2E.  ♦  A 
3y  5 


32F 

3y2 
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in  which  F,  where  T  -~ 

F*  =  {U,  W,  T,  q^} , 


(18) 


(19) 


is  the  vector  of  the  principal  dependent  variables,  A^  through  A^  are  square 
matrices  of  order  4,  and  A^  is  a  four-dimensional  vector.  A^  through  A^  are 
functions  F,  but  are  regarded  as  known,  at  any  given  iteration  level,  via 
the  customary  linerization  process. 


Equation  (18)  is  solved  in  a  three-dimensional  domain:  two  space  dimensions, 
x,y  describing  a  plane  normal  to  the  generators  of  the  cylinder,  and  the 
time  dimension,  t.  The  calculation  advances  in  the  positive  time  direction 
yielding  the  complete  solution  in  the  x  -  y  plane  at  each  step  (Figure  1). 

The  node  points  in  the  x  -  y  plane  are  scanned  by  an  alternating-direction 
(ADI)  technique:  a  partial  solution  is  first  sought  along  lines  x  =  constant, 
and  a  revised  solution  is  then  sought  along  lines  y  =  constant,  (Figure  2). 

One  complete  iteration  consists  of  a  scan  in  each  direction. 

There  are  several  advantages  of  the  ADI  technique:  the  number  of  iterations 
required  is  a  relatively  weak  function  of  the  node-point  density,  the 
stability  of  the  scheme  is  not  compromised  by  ambiguities  in  sign  of  the 
chordwise  velocity  component:  U,  and  extension  to  higher-order  systems  of 
equations  is  readily  accomplished  (an  ADI  scheme  is  used  in  the  method 
developed  at  Sybucon  for  integrating  the  unsteady  Navier-Stokes  equations; 
see  Reference  9) . 
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For  a  scan  in  either  direction.  Equation  (18)  takes  the  form 

PG  =  Q  (20) 

where  the  G  is  the  vector  whose  elements  are  the  values  of  F  at  each  node 
point  along  the  line  of  scan.  P  is  a  tridiagonal  matrix  whose  elements  are 
4x4  submatrices,  and  Q  is  a  vector  whose  elements  are  themselves  4x1 
matrices . 

Second-order  differencing  is  used  throughout  the  method.  Backward  differences 
(relative  to  the  local  flow  direction)  are  used  for  the  diagonal  elements  of 
and  A^,  in  Equation  (18),  and  central  differences  for  the  diagonal  elements 
of  Aj_  and  the  off-diagonal  elements  of  A^  (A^  and  A,.  have  no  non-zero  off- 
diagonal  elements).  In  order  to  preserve  the  second-order  status  of  the 
backward  differences,  two  node  points,  at  the  current  time  level,  and  three 
at  the  previous  time  level,  are  involved  in  each  representation  of  a  derivative. 
The  differencing  arrangement  is  chosen  primarily  for  stability  reasons,  but  it 
offers  good  precision  characteristics  also.  The  solution  of  Equation  (20): 

G  +  P_1Q  (21) 

yields  the  solution  for  the  particular  scan  and  iteration. 
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Ill  THE  VALIDATION  OF  THE  METHOD 


I 


A  range  of  calculations  was  performed  to  establish  the  credibility  of 
the  method.  For  two-dimensional  laminar  flow  a  time -relaxation  solution 
was  found  to  be  in  excellent  agreement  with  the  classical  Blasius  solution 
(Figure  3) .  Numerical  solutions  were  derived  for  various  node-point 
densities,  as  shown,  and  it  is  evident  that  satisfactory  convergence  is 
being  achieved  with  second-order  accuracy  (Figure  4);  an  error  in  wall 
shear  stress  of  less  than  0.2%  is  incurred  with  as  few  as  11  vertical 
node-points.  Some  comparison  calculations  show  that  the  corresponding 
first-order  computations  result  in  substantially  larger  errors  (Figures  3,41. 

A  time-relaxation  calculation  for  a  steady  turbulent  flow  is  shown  in 
Figures  5,6,  and  a  typical  level  of  accuracy  is  achieved  with  respect  to 
experimental  data.  With  regard  to  the  validation  of  the  turbulence  model 
is  is  worth  mentioning  that  a  separate  study  has  been  conducted  to  verify 
a  three-dimensional  "steady"  version  of  the  present  method  [13].  The  results 
of  this  study  showed  that  some  small  improvements  could  be  made  to  the 
empirical  content  of  the  model  although,  in  general,  the  model  performed 
reasonably  well.  The  corresponding  improvements  were  not  included  in  the 
present  method  because  of  doubts  about  the  invariance  of  the  proposed 
changes  to  rotation  of  the  coordinate  system. 


In  order  to  examine  the  effect  of  grid  density  on  the  solution,  under 
reversed-flow  conditions,  some  calculations  were  done  for  a  retarded 
turbulent  boundary  layer,  and  the  results  are  presented  in  Figures  7  through 


10.  The  flow  is  similar  to  those  considered  in  Sections  III  through  V, 
below.  The  comparisons  show  that  the  details  of  the  flow  are  represented 
surprisingly  well  even  with  only  11  vertical  node-points.  The  remaining 
calculations,  presented  throughout  this  report,  were  done  using  21  points 
in  the  vertical  direction. 
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IV  CALCULATIONS  FOR  TWO-DIMENSIONAL  INCOMPRESSIBLE  FLOW 


The  calculations  were  all  performed  for  the  class  of  prescribed  external 

g 

velocity  distributions  described  as  "frozen  flows"  .  Specifically 

Ug  =  Uq,  for  t  <_  o  and  all  x,  (22) 

U 

jp  =  1  -  —  {1-  f(t)},  for  t  >  o  and  o  <  x  £  x. ,  (23) 

o  1  1 

U  x  -x 

7j—  =  1  -  -  {1  -  f(t)},  for  t  ^  o  and  x.  <_  x  <_  c.  (24) 

o  X1  1 

The  function  f(t),  and  the  form  of  the  velocity  distribution,  are  shown  in 
Figure  11.  This  external  velocity  distribution  imposes  an  adverse  pressure 
gradient  (whose  severity  increases  with  time)  over  the  forward  part  of  the 
chord,  and  a  favorable  gradient  over  the  rearward  part. 

Figures  12,  13  show  typical  solutions,  for  this  type  of  velocity  field,  for 
both  laminar  and  turbulent  incompressible  flow.  It  will  be  seen  that  the 
passage  through  the  point  of  zero  wall  shear  stress  is  smooth  and  uneventful, 
and  that  substantial  regions  of  reversed  flow  develop  before  the  onset  of 
singular  conditions  occurs.  Figure  14  shows  the  rapid  thickening  of  the 
boundary  layer  which  takes  place  some  time  after  the  first  appearance  of 
reversal.  For  the  particular  combination  of  parameters,  in  Equations  (22-24), 
singularity  onset  occurred  before  the  freezing  of  the  external  flow.  Other 
calculations,  in  which  these  parameters  were  different,  showed  the  approach 
to  singular  conditions  after  the  external  flow  was  frozen,  but  the  diver¬ 
gence  of  the  solution  was  essentially  the  same.  Figures  15,  16  show  the 


familiar  spacial  distributions  of  displacement  thickness  and  wall  shear 
stress,  for  incompressible  flow.  These  are  similar  to  the  distributions 
presented  in  several  earlier  papers  [5,6,9]. 

Figure  17  shows  the  effect  on  the  solution  of  allowing  the  dissipation 
length  to  lag  behind  the  local  boundary- layer  thickness.  It  was  noted 
earlier  that  the  inclusion  of  the  rate  equation  for  dissipation  lenght 
raises  the  turbulence  model  to  two-equation  status.  As  would  be  expected, 
the  effect  of  the  lag  is  most  conspicuous  in  regions  of  rapid  boundary- 
layer  thickening.  On  the  other  hand,  the  influence  on  the  velocity  and 
turbulent-kinetic-energy  fields  is  relatively  small.  Figure  18  shows  the 
results  of  making  a  rough  correction  for  the  effects  of  mementum  transport 
normal  to  the  surface.  Here  the  normal  pressure  gradients  were  calculated 
retrospectively  after  each  iteration,  and  were  averaged,  across  the  boundary 
layer.  In  subsequent  iterations  of  the  solution  the  effect  was  reintro¬ 
duced  in  the  form  of  a  y-dependent  modification  of  the  streamwise  gra¬ 
dients.  The  procedure  appeared  to  be  stable,  and  resulted  in  a  small 
change  in  the  converged  solution. 
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V  CALCULATIONS  FOR  YAWED  FLOW  AND  FOR  COMPRESSIBLE  FLOW 


The  effect  of  yaw  on  the  solution  was  investigated  by  adding  a  constant 
W-component  to  the  prescribed  external  velocity  field.  The  results 
presented  in  Figures  19,  20,  which  are  typical  of  those  obtained,  show 
that  the  effect  on  the  chordwise  velocity  field  within  the  boundary  layer 
is  relatively  minor.  The  spanwise  component  of  wall  shear  stress  does 
not,  of  course,  reverse  although  it  exhibits  a  shallow  minimum  near  the 
first  point  of  reversal  of  the  chordwise  component.  The  angle  between 
the  external  streamlines  and  the  x-axis  varies  with  x,  reaching  a  maxi¬ 
mum  at  x  =  Xj.  At  the  time  level  represented  in  Figure  20,  this  maximum 

angle  is  63°,  while  the  maximum  angle  of  the  wall  streamlines  is  in  the 
neighborhood  of  135°. 

Figures  21,  22,  show  the  effect  of  compressibility  on  the  solution.  It 
will  be  seen  that  there  is  a  striking  change  in  the  distribution  of  dis¬ 
placement  thickness  as  the  Mach  number  increases.  At  a  Mach  number  of  3 
there  is  little  sign  of  the  rapid  boundary- layer  thickening  which  domi¬ 
nates  the  picture  at  low  speeds.  Another,  less  conspicuous,  effect  of 
increasing  Mach  number  is  the  delay  of  reversal  onset;  this  is  the  only 
significant  effect  on  the  distribution  of  wall  shear  stress  in  the  neighbor¬ 
hood  of  reversal.  Interestingly,  though,  these  calculations  indicate 
little  variation  of  the  position  of  reattachment. 
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VI  OTHER  CALCULATIONS 


A  number  of  important  calculations  were  performed  to  investigate  the 
effects  of  alleviating  the  external  velocity  gradients  so  as  to  delay 
the  onset  of  singular  conditions.  These  calculations  have  been  reported 
elsewhere  [12],  and  will  not  be  presented  in  detail  herein. 

Briefly,  it  was  demonstrated  that  by  judicious  alleviation  of  the  chord- 

wise  gradients  the  rapid  boundary- layer  thickening,  indicative  of  incipient 

singular  conditions,  could  be  delayed  indefinitely.  The  intention  was  to 

simulate  the  modification  of  the  pressure  field  which  occurs  in  a  real  flow 

as  the  result  of  coupling  between  the  boundary  layer  and  the  outer  flow. 

Typical  results  are  reviewed  in  Figures  23-28,  for  laminar  as  well  as  for 

turbulent  flow.  It  will  be  seen  that  the  alleviation  of  the  gradients  is 

confined  to  the  regions  where  the  singularity  would  otherwise  be  expected 
to  develop.  With  the  alleviation  mechanism  in  operation,  the  calculations 

could  be  continued  far  beyond  the  point  where  the  calculations  in  pre¬ 
scribed  external  velocity  fields  broke  down.  The  rapid  boundary- layer 
thickening  was  contained,  and  the  region  of  reversed  flow  remained  of 
roughly  constant  extent  as  the  solution  asymptoted  towards  steady-state 
conditions. 

As  was  noted  in  Reference  12,  these  calculations  showed  that  the  singularity 
onset  is  not  an  inevitable  feature  of  first-order  boundary- layer  theory,  but 
appears  to  result  from  an  inappropriate  use  of  the  theory:  i.e.,  in  pre¬ 
scribed  external  gradients  which  do  not  properly  allow  for  the  effects  of 
interaction  between  the  visious  and  inviscid  regions. 
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VII  CONCLUDING  REMARKS 


The  purpose  of  this  research  was  two-fold:  to  develop  a  more  compre¬ 
hensive  computation  tool  for  the  particular  class  of  flows  than  was 
hitherto  available,  and  to  demonstrate  the  use  of  the  method  as  part  of  an 
ongoing  study  of  specific  fluid-dynamic  phenomena. 

The  new  calculation  method  appears  to  be  firmly  based,  both  from  the  stand¬ 
point  of  the  underlying  physical  model  and  in  terms  of  the  numerical  tech¬ 
niques  used  to  integrate  the  governing  equations.  Although  the  validation 
of  the  method  has  not  been  extensive,  tests  run  to  date  indicate  that 
acceptable  accuracy  and  reliability  can  be  achieved  at  reasonable  compu¬ 
tation  cost. 

The  fluid-dynamic  studies  conducted  by  the  new  calculation  method  included 
parallel  computations  of  unsteady  laminar  and  turbulent  flows  with  embedded 
regions  of  reversal.  A  broad  measure  of  similarity  emerged  from  the  solu¬ 
tions  in  the  two  types  of  flow.  As  was  to  be  expected  from  the  results  of 
earlier  work,  the  development  of  the  solution  for  either  type  of  flow,  was 
interrupted  by  the  onset  of  a  singularity  some  time  after  reversal  was 
observed. 

Some  calculations  were  performed  to  investigate  the  effects  of  Mach  number, 
and  the  effects  of  yaw,  on  the  development  of  reversed  flows.  Yaw  seemed 
to  have  relatively  little  effect  on  the  development  of  the  chordwise  flow, 
but,  of  course,  gave  rise  to  spanwise  flow  components  which  have  their  own 
interesting  dynamics.  The  principal  effect  of  increasing  Mach  number  was 
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o  reduce  the  thickness  of  the  boundary  layer  in  the  neighborhood  of  reversal; 
at  supersonic  speeds  only  a  limited  build-up  of  displacement  thickness  was 
observed,  whereas  a  rapid  build-up  is  familiar  in  incompressible  flow. 


A  significant  outcome  of  the  present  research  has  been  the  observation  that 
onset  of  the  separation  singularity  in  unsteady  flow  can  be  delayed,  or 
avoided  altogether,  by  alleviation  of  the  gradients  in  the  external  flow. 
The  alleviation  of  the  gradients  was  intended  to  simulate  the  interaction 
between  the  boundary  layer  and  the  outer  flow.  The  calculations  indicated 
that  the  singularity  onset  could  be  delayed  more-or-less  indefinitely,  and 
that  conditions  of  steady-state  reversal  could  be  approached  by  way  of  an 
asymptotic  calculation  for  long  times.* 

The  new  calculation  methoc  will  presumably  be  used  for  continued  basic 
research  on  the  mechanics  of  unsteady  reversed  flows.  In  addition  the 
method  is  currently  being  used  as  a  component  tool  in  a  study  of  unsteady 
flow  over  finite  wings.  This  latter  work  is  part  of  an  on-going  program 
being  supported  by  the  Office  of  Naval  Research. 


*  Here,  "long"  refers  to  the  physical  time  represented  rather 
than  to  the  required  computation  time. 
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partial  implicit  solution 
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Figure  2  SCANNING  OF  THE  x  -  y  PLANE  USING  THE  ALTERNATING- 

DIRECTION  TECHNIQUE 
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Figure  3  TIME-RELAXATION  SOLUTION:  2-DIMENSIONAL  LAMINAR  FLOW 
OVER  A  FLAT  PLATE 
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Figure  5  TIME  RELAXATION  SOLUTION:  2-DIMENSIONAL  TURBULENT  FLOW 
OVER  A  FLAT  PLATE 


Figure  7  EFFECT  OF  VARIATION  OF  VERTICAL  GRID  DENSITY 
[UQt/c  =1.3,  x/c  =  0.66] 


Figure  8  EFFECT  OF  VERTICAL  GRID  DENSITY 
[UQt/c  ■  1.3,  x/c  *  0.66] 
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Figure  12  LAMINAR  BOUNDARY  LAYER  PENETRATING  A  REVERSED-  FLOW  REGION 
(Laminar  Frozen  Flow:  Rg  =  104,  =  0.5,  =  2.0  ,  x/c  =  0.5) 
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Figure  13  TURBULENT  BOUNDARY-LAYER  CALCULATION  PENETRATING 
A  REVERSED-FLOW  REGION  (Turbulent  Frozen  Flow:  Rg  =  107,  ff  =  0.25 

=  1.5,  x/c  =  0.66) 
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Figure  14  RAPID  THICKENING  OF  THE  BOUNDARY  LAYER  SUBSEQUENT 

TO  REVERSAL  ONSET 
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Figure  16  SPACIAL  DISTRIBUTION  OF  WALL  SHEAR  STRESS 
[M  =  0,  2-Dimensional  Turbulent  Flow] 
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Figure  17  COMPARISON  BETWEEN  1-  AND  2-EQUATION  TURBULENCE  MODEL 
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1st  order  theory 


Figure  20  EFFECT  OF  YAW  ON  WALL  SHEAR  STRESS 
[M  =  0,  Turbulent  Flow] 


EFFECT  OF  MACH  NUMBER  ON  DISPLACEMENT  THICKNESS 
[2-Dimensional  Turbulent  Flow] 


[2-Dimensional  Turbulent  Flow] 


U  t/c 


Figure  23 


Figure  24  LAMINAR  FLOW  KITH  ALLEVIATION  OF  PRESSURE  GRADIENTS 

DISPLACEMENT  THICKNESS 
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OF  ALLEVIATION:  IROZEN  TURBULENT  FLOW 


Figure  21  TURBULENT  FLOW  WITH  ALLEVIATION  OF  PRESSURE  GRADIENTS: 

DISPLACEMENT  THICKNESS 


TURBULENT  ELOW  WITH  ALLEVIATION  OF  PRESSURE  GRADIENTS 
WALL  SHEAR  STRESS 


